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Aims. Recent theoretical predictions for the winds of Wolf-Rayet stars indicate that their mass-loss rates scale with the initial 
stellar metallicity in the local Universe. We aim to investigate how this predicted dependence affects the models of Wolf-Rayet 
stars and their progeny in different chemical environments. 

Methods. We compute models of stellar structure and evolution for Wolf-Rayet stars for different initial metallicities, and 
investigate how the scaling of the Wolf-Rayet mass-loss rates affects the final masses, the lifetimes of the WN and WC subtypes, 
and how the ratio of the two populations vary with metallicity. 

Results. We find significant effects of metallicity dependent mass-loss rates for Wolf-Rayet stars. For models that include 
the scaling of the mass-loss rate with initial metallicity, all WR stars become neutron stars rather than black holes at twice 
the solar metallicity; at lower Z, black holes have larger masses. We also show that our models that include the mass-loss 
metallicity scaling closely reproduce the observed decrease of the relative population of WC over WN stars at low metallicities. 
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Massive stars are dominant sources of energy and nucle- 
osynthesis products for the interstellar medium (ISM). 
This input is provided during their entire lives via stel- 
lar winds, followed by spectacular deaths as supernovae 
and long-duration gamma-ray bursts (GRBs). Because of 
their intrinsic brightness, they can be identified individu- 
ally in extragalactic galaxies, whilst their integrated emis- 
sion dominates the spectra of distant galaxies in the high 
redshift Universe. Wolf-Rayet (WR) stars are bare helium 
stars, with spectra that are characterised by broad emis- 
sion lines, and a lack of hydrogen (H), which is a result of 
mass loss. Stellar winds continue to enrich the ISM dur- 
ing the WR phase, where the objects turn into important 
sources of helium (He), nitrogen (N), carbon (C), oxygen 
(O), and other elements in our Universe. In addition, the 
N-rich WN and C-rich WC stars are the favoured pro- 

J :enitors of type Ibc supe rnovae due to their lack of H 
Ensman fc Wooslevlll988(l : more recently they have be- 
come suspect to be th e progenitors of long-duration GRBs 
l)Popham et al.lll999|) . 

Nonetheless, WR stars remain somewhat of an enigma. 
Measuring the exact mass-loss rates (M) and other stel- 
lar parameters is particularly challenging, as WR atmo- 
spheres are not in local thermodynamic nor hydrostatic 



equilibrium, which prevents a direct comparison of surface 
temperatures and radii with stellar models. WR winds are 
optically thick and inhomogeneous, and wind-clumping 
has led to severe overestima tions of their mass-loss rates 
by factors of about t hree ijHamann fc Koesterkd llQQSt 
iNueis fc Lamersll2000ll . The intricacies of clumping and 
their effect on measurements of WR mass-loss rates re- 
main a source of uncertainty today. 

Despite these uncertainties, we have significantly ad- 
vanced our understanding of WR stars, as it has be- 
come increasingly c lear that W R winds are driven by 
radiation pressure Lucv fc Abb ott 1993; Hillicr 2 003t 
iGrafener fc HamannI |2005|) . and that WR mass loss 
dep ends on their initial metallicity {Z). In particu- 
lar, ICrowther et al. ( 2002J presented observational evi- 
dence for a WR M — Z relation for WC stars, while 
IVink fc de 12005") recently showed how M is pre- 

dicted to vary with initial metallicity for late-type WC as 
well as WN stars. 

In this article, we investigate the implications of this 
AI — Z scaling for WN and WC stars. First we de- 
scribe the construction of our stellar models and the 
mass- loss scheme we have implemented (Sect. EJ. In 
Sect, m we discuss how o ur models compare to models of 
iMevnet fc Maede^ l)2005(l . We subsequently predict final 
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Fig. 1. Hertzsprung- Russell diagram comparison be- 
tween our solar metall icity models and those from 
iMevnet k, Maedeil l)2003l) . Our models are represented 
by the s olid grey lines, the other lines are taken from 
iMevnet fc Maedcr ( 2003). The dashed lines indicate non- 
rotating models, while the dotted lines represent an initial 
rotational velocity of 300 kms~^. The stars in the upper 
panel have an initial mass of 25AfQ, while those in the 
lower panel have an initial mass of AQMq. 



masses (Sect.^ and WR lifetimes (Sect.E)), and compare 
the model results to the observed WC/WN ratio at vari- 
ous metallicities (Sect.|Sl). We finally check the sensitivity 
of our results to potential remaining uncertainties in the 
exponent of the scaling of WR mass-loss rates, before we 
conclude in Sect.H 



2. Construction the Stellar Models 

Our stellar models are produced with the Cambridge 
stellar evolutio n code, STARS, originally de veloped by 
lEggletoiil lll97lh and upda ted most recently bv lPols et al.1 
l|l995|) and lEldridge &: Tbut (,2004a|l . Further details can 
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Fig. 2. Same as Fig.^for initial masses of 60AfQ (upper 
panel) and I2OM0 (lower panel). 



be found at the code's home page ^ . The models are similar 
to the o vershooting models described in Eldridge & Tout! 
((2OO41J 1. The differences are as follows: the initial masses 
are 25, 30, 35, 40, 45, 50, 55, 60, 70, 80, 100 and 120Mq. Z, 
the initial metallicity, takes the values 0.001, 0.002, 0.003, 
0.004, 0.008, 0.02 and 0.04, i.e. the range is ^Zq to 2Zq. 
The metallicity comprises scaled-solar abundances. In con- 
structing the stellar models, we assume "solar metallicity" 
to correspond to a metal mass fraction Z ~ 0.02 and scale 
the mass-loss rates from this point. The most r ecent anal- 
ysis of the Sun's composition (Asolund et af] l2005j) sug- 
gests that this metallicity may need to be revised to a 
value as low as Z = 0.012. Of course the metallicity of 
the solar neighbourhood may not exactly equal Zq itself, 
and it is the metallicity of the solar neighbourhood for 
which mass-loss rates have been measured and calibrated 
against. Furthermore, the relative scale is more important 
than the absolute value, and the uncertainties related to 



^ http://www.ast.cam.ac.uk/~stars 
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the value of the initial metallicity are well within the un- 
certainty of WR mass-loss rates. 

The mass-loss scheme is the most relevant factor in 
constructing the stellar models. For OB s tars, we use the 
mass-loss predictions of I Vink et al.l l)200l|) ^ which include 
a metallicity sca ling. For all ot her pre-WR phases, we em- 
ploy the rates o fEZ Jager et al. (1988) scaled with metal- 
licity by a factor of (Z/Zq)^-^. When the star becomes 
a WR star (X^.^rf^r^ < 0.4, l ogfTpff /JC) > 4.0), we use 
the rates of iNugis fc Lamera ()2000|) . It is the scaling of 
these rates with Z that produces different sets of mod- 
els. The scaling we apply depends both on the WR star 
type and the initial metallicity of the star. The WR star 
type is determined as follows. The star is initially a WNL 
star, and once Xgurfaco = 0, the star becomes a WNE star. 
To determine the switch from WNE to WC star, we de- 
fine the following parameter, ( = {xq + xo)/y, where xc, 
xo and y are the surface number fractions of C, O and 
He. Following Macder & Mevnet (1994), the star becomes 
WC when Xgurface = and ^ > 0.03. In a similar vein, 
the switch to a WO star occurs when C > 1.0, but we 
only encounter the WO phase when we do not include the 
AI — Z scaling. Since WO lifetimes are much shorter than 
those of the other WR phases, we do not distinguish WO 
stars from WC stars. 

With the WR type defined, we apply the M scal- 
ing factor, (Z/Zq)'^, where x is shown in Tabled When 
Z < 0.002 the scaling factor changes to (Z/0.002)^2 x 
(0.002/^0)^1 achieve the correct scaling, xi is the ex- 
ponent when Z > 0.002 and X2 is the exponent when 
Z < 0.002. In scheme A, M does not scale with Z. The 
preferred values of our M — Z dependence (scheme B) 
are taken fro m the study of late-type WN and WC M 
predictions of IVink fc de Koter (2005) , where the expo- 
nents come from those predicted for the wind momentum 
(M Voo) dependence on Z. In other words, the underly- 
ing assumption for scheme B is that the wind velocity of 
WR stars does not vary with Z , as we anticipate a po- 
tential Z-dependence of Voa to be relatively weak. The 
reason is that the line force due to iron lines is predom- 
inantly found to deter mine the mass-loss rate, whilst th e 
CNO elements set v^o ()Vink et al.lll999HPuls et al.ll200(]() . 
Because these CNO abundances are largely determined by 
"self-enrichment" , the Voo versus initial stellar metallicity 
dependence (which scales with Fe) is expected to be even 
weaker for WR stars than it is for OB stars. Nonetheless, 
to accommodate for a potential metallicity dependence of 
the WR wind velocity, we include scheme C, where the 
exponents of scheme B have been slightly reduced. This 
scheme will be used for comparison later as to estimate 
the importance of the potential wind velocity metallicity 
dependence on WR models. We finally include scheme D 
to relate our computed exponents to the commonly quoted 
constant exponent of 0.5. 



Table 1. Scaling of the WR mass-loss rates with metal- 
licity. The values represent the exponent of M{Z) = 





z 


Scheme A 


B 


C 


D 


WN 




0.0 


0.85 


0.69 


0.5 


WC 


>0.02 


0.0 


0.4 


0.3 


0.5 




<0.02 


0.0 


0.66 


0.56 


0.5 




<0.002 


0.0 


0.35 


0.25 


0.5 



^ An IDL routine to compute the mass-loss rate as 
a function of stellar parameters is publicly available at 
http : //www . astro . keele . ac . uk/ ~ j sv/ 



3. Model comparison 

We compare our stellar models to other contemporary 
models of WR stars. We do this for four initial masses in 
Figs, n and [21 using the n on-rotating and rotating mod- 
els of Mevne t fc Maeded ({2003) and iMevnet fc Maederl 
(I2OO5I) , the Geneva models. The surface temperatures have 
not been corrected for the no n-equilibrium effects of th e 
WR atmosphere as described in lMevnet fc Maedeil 1 20051) . 
Although stellar rotation is inherently a 3D process, the 
Geneva models include the effects of stellar rotation by 
using approximations that can be included within a ID 
evolution code. We note that the Geneva tracks end after 
core helium burning, while ours progress onto core carbon 
burning. 

Comparing the 25M0 models in Fig.^ we find that the 
non-rotating Geneva model agrees well with our model, al- 
though the Geneva track does not include the RSG phase. 
The rotating model differs significantly in that the star 
ends its evolution as a WR star. This is due to a combina- 
tion of rotation enhancing the mass-loss rates, rotation in- 
creasing the main-sequence lifetime leading to extra mass- 
loss, and a larger helium core due to the rotational mix- 
ing. The 4OM0 models all end as WR stars. However, the 
Geneva tracks never become cool red supergiants as ours, 
and the rotating model produces more luminous stars due 
to enhanced energy production as a result of rotational 
mixing. The final masses of both types of Geneva models 
are a few Mq more massive than ours. Similar results are 
found in Fig. [5]for the 6OM0 models, but in addition the 
rotating model never moves over towards the red side of 
the HR-diagram. For the I2OM0 models, stars all remain 
in the blue part of the HR-diagram for their entire evo- 
lution. The rotating star in this cases achieves the same 
masses as in our model. 

Overall the differences are relatively minor. They can 
be explained by tw o important differences b etween our 
models and those of iMevnet fc Maedeil ()2003(l . Most de- 
tails of the evolution codes are however very similar. For 
example, both include a small amount of convective over- 
shooting. While both rnqdels employ the same base m ass- 
loss rates oflVink et al.l ()200l|) . ldeJager et all (|l988l) and 
iNugis fc Lamerj (|2000|) : their use is different. In our mod- 
els we use all of the rates unm odified as d e tailed above. 
The Geneva models using the IVink et all l|200lh rates, 
which are for non-rotating stars as a basis to reduce the 
empirical rates to account for the rotation of the observed 
stars. They then enhance the mass-loss rate according to 
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the rotation rate of the model. This leads to the mass-loss 
rate for their non-rotating models being less than those of 
our models, whilst the mass-loss rates employed for their 
rotating models are very similar to those of our models, 
as we use unaltered mass-loss rates. The above implies 
that our models are similar to the rotating Geneva mod- 
els, except that we do not consider the effects of rota- 
tional mixing. The extra mixing of hydrogen into the stel- 
lar core extends the main-sequence lifetime for rotating 
stars, which leads to smaller post main-sequence stellar 
masses. In addition, mixing leads to greater helium core 
masses, so that larger WR star masses are obtained in the 
Geneva models than in ours. The Geneva non-rotation 
models have a higher minimum initial mass for WR stars, 
37 Mq, compared to our 27 Mq, and shorter WR lifetimes. 
While the rotating models have a lower minimum initial 
mass of 22Mq and longer lifetimes. It is difficult to make a 
comparison at other metallicities as there are no contem- 
porary non-rotating Geneva models at lower metallicities. 
Nonetheless, similar trends are found comparing the ro- 
tating models to our scheme A models. While the scheme 
B models differ by producing much greater final masses 
than do the same Geneva models. 

In summary, the models we present are non-rotating 
models that have not been modified i n anyw ay to account 
for rotation as in iMevnet fc Maede^ l)2003|) . Importantly 
we do not consider how rotation introduces enhanced mix- 
ing and we do not alter the mass-loss rates to account for 
the rotation of the stars. 

4. Final WR star masses 

We first consider how relevant aspects of the stellar mod- 
els are affected by the scaling of the WR mass-loss predic- 
tions. The final mass of our stars is directly determined by 
the amount of mass lost over their lifetimes. We present 
our final mass predictions in Fig. |21 where we show results 
at four metallicities: solar, twice solar, two fifths solar and 
one fifth solar. In the upper panel for scheme A we see 
that lower metallicity WR stars are more massive due to 
the reduced pre-WR mass loss so the initial WR masses 
are larger. In the solar case, stars with initial masses > 
27Mq become WR stars. As the initial mass increases, so 
does the final mass up to a peak at 8OM0, after which, the 
final mass again decreases. This is due to the time spent 
as a red supergiant becoming shorter, and effectively, the 
stars remain hot for their entire lifetimes, experiencing 
strong mass-loss during their main-sequence lives leading 
to lower masses. 

Comparing these results to the lower panel of Fig. 13 
we find large differences. At twice solar metallicity the 
masses are decreased by a few solar masses. While at lower 
metallicities the final masses are increased by the same 
factor two fifths solar. Furthermore, the change in final 
mass from solar to these LMC-type mode ls are similar to 
those predicted bv ICrowther et all l(2002l) . At the lowest 
metallicity here the effects are greatest and we find the 
masses to be larger by around lOM©. 
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Fig. 3. The final masses from our models versus the ini- 
tial masses at four different metallicities. The upper panel 
shows the predictions according to scheme A, where no 
WR scaling has been applied. The lower panel denotes 
predictions according to scheme B where the WR scaling 
is applied. 

The effects on the final masses of WR stars are most 
relevant for the prediction of the masses of the compact 
remnants. If we assume that stars with He cores greater 
than %Mpi p roduce black holes at the end of their lives 
l)Heger et alJl2003J . we find that by including the WR 
scaling at twice solar metallicity, all compact remnants 
are predicted to be neutron stars rather than black holes. 
On the other hand, at lower Z, the masses of black holes 
are predicted to be much larger with the scaling. In other 
words, the compact remnant mass function shows a strong 
Z-dependence. 

5. WR star lifetimes 

The lifetimes of our WR models depend on the mass-loss 
rate in two ways. First, M determines the total lifespan of 
the WR star by affecting the total stellar mass and thus 
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Fig. 4. The variation of the WR hfetimes with initial 
mass when the mass- loss rates are not scaled (scheme A). 
Each region is labelled but progresses from bottom to top 
by WNL, WNE and WC phase. The upper panel is for 
Z=0.04 and the lower panel is for Z=0.02. 



the vigour of the nuclear reactions in the core. Second, 
stronger mass loss strips mass quicker off the stellar sur- 
face, leading to the exposure of the products of later nu- 
clear reactions, resulting in the appearance of different 
WR subtypes. 

Figures ^ and |S1 show the times spent in the various 
phases when we do not scale M with Z (Scheme A). At 
higher metallicities, the WR lifetime is longer, because of 
larger pre-WR mass loss so that the WR phase is entered 
at an earlier stage resulting in a longer period of core he- 
lium burning during the WR phase. The kinks seen in 
the solar and twice solar Z plots (Figs. 0] and are due 
the fact that stars with masses below the kink undergo a 
red supergiant phase, where mass is stripped and the sur- 
face H abundance drops below the limit set for WNL stars 
(Xsurfacc = 0.4). Howcvcr, the stars are still cool and some 
time passes before the stars are hotter than the tempera- 
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Fig. 5. Same as Fig. EJ but with metallicities Z=0.008 
(upper panel) and Z=0.004 (lower panel). 



ture limit (10^ K). For stars more massive than the kink, 
the stars spend most of their evolution at temperatures 
hotter than this limit and become WR stars as soon as 
the surface abundance requirement is reached. 

Figures and [71 show the times spent in the various 
WR phases when the predicted mass-loss scaling is in- 
cluded (Scheme B). The differences between the Z-scaled 
and unsealed models are relatively minor as far as the 
total WR lifetimes are concerned, although bX Z ~ 0.04 
the lifetimes are somewhat longer in scheme B, as the WR 
stars are less massive and the nuclear reactions proceed at 
a slower rate, extending the lifetime. The reverse occurs 
at sub-solar Z . 

When we consider the lifetimes of the WR sub-phases, 
the mass-loss scaling affects how quickly nuclear burning 
products are exposed at the surface. At Z = 0.04, we find 
that for stars below SOMq the WN lifetime is reduced, 
while it is increased for stars above this value. The en- 
hanced mass-loss leads to smaller He burning cores, so a 
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Fig. 6. Same as Fig.01but with the mass-loss rates scaled 
according to scheme B. 

longer time must pass before the He burning products are 
exposed at the stellar surface. 

At Z = 0.008 and Z = 0.004 wc find that the lower 
mass objects end their lives as WNL stars and do not 
progress any further along the WR evolutionary path. 
This is due to the weakening mass-loss stripping of a large 
fraction, but not all of the H. For more massive stars it 
leads to a reduction of the WC lifetime to the benefit of 
the WN lifetime. At even lower metallicities, the M — Z 
scaling leads to a complete absence of WC stars in our 
single star computations. 

6. WC to WN populations ratio versus metallicity 

To evaluate whether the predicted M — Z scaling is simi- 
lar to that in nature, we wish to compare our predictions 
to observations. We have compared the observed mini- 
mum initial mass r equired for the different WR phases 
l|Massev et al.ll200(]|) as well as the o bserved ratio of type 
Ibcto type II SNe versus metallicity l|Prantzos fc BoissieJ 
l2003l) to our predicted values. We find that the inclu- 
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Fig. 7. Same as Fig.Elbut with the mass-loss rates scaled 
according to scheme B. 



sion of the predicted mass-loss scaling provides a better 
agreement with these observations, however the useful- 
ness of these tests is limited by their large uncertainties. 
Fortunately, a more robust test is possible: the ratio of the 
WC versus WN star populations at different metallicities. 
We weig ht the WR lifetimes (Sect. |S} by an initial mass 
function l|KrouDa et al...l99Sl and predict WC/WN ratios 
that we compare to observed values. We implicitly assume 
a constant star-formation rate. 

FigurelHlshows these predictions as well as the observed 
values. For the observed values the metallicity mass frac- 
tion was calculated from the log]^o(0/H) -I- 12 values by 
comparing them to the values from our sets of models. 
Because this process is ambiguous and as discussed above 
the position of solar metallicity for the mass-loss scaling is 
indistinct, we ha ve assumed the m etallicities are uncertain 
by 25 percent. In lMevnet fc Maed er (200^ a similar graph 
is presented where the x-axis was left as log]^Q(0/H) -I- 12. 
In that plot an assumed solar value was taken rather than 
the composition of the stellar models. We here re-plot the 
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Fig. 8. The WC/WN ratio versus me t allici t 
Observed po ints are fr om iMassev fc JohnsoiJ 1199^, . 
iMassev fc Duf fy ( 2 00l|l : iMassev fc Holmes! ((2002j) . Tlie 
error bars for the metallicity of the observed points reflect 
the uncertainty of converting from logig(0/H) + 12 and 
the uncertainty in the exact value of solar metallicity and 
therefore where the mass-loss rates should be scaled from. 
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Fig. 9. Same as Fig. |Hlbut using the Geneva models. 



Geneva model results in Fig. |51 on the same scale as in 
Fig. El 

Despite the uncertainties in the observations it is clear 
that the obser ved WC/WN ratio decre ases at lower Z. 
The models of iMevnet fc Maedej l)2005|) suggest that by 
mixing a population of rotating and non-rotating stars, 
agreement can be made between their models and the 
observations (Fig. O. In Fig. |Slwe show that agreement 
can be reached without any such assumption. All that is 
required is to scale the WR mass-loss rates with initial 
metallicity. 

A complete evaluation of the relative importance of 
our WR mass-loss scaled models versus models including 
stellar rotation is difficult, as there are no contemporary 
non-rotating Geneva models available at lower metallici- 
ties. In addition, the Geneva models only include a WR 
metallicity scaling for supersolar models, despite the fact 
that the evidence used for this supersolar scaling comes 



from observations of WR stars in the low metallicity en- 
vironment of the LMC. We nevertheless anticipate that 
including a WR M — Z scaling in rotating models will lead 
to an increase in the final WR masses (predicted by the 
Geneva code) and a decrease of the WC lifetimes, possi- 
bly even lowering the WC/WN ratio to values below those 
observed. 

In other words, Fig.|Slhints that as far as the WC/WN 
ratio is concerned, rotation is of secondary importance to 
employing the correct WR scaling of the mass-loss rates 
with metallicity. This assertion does not imply that rota- 
tion is not important for massive star evolution, as there 
is significant observational evidence for rotation-induced 
mixing (e.g. Howarth & Smith 2001), but we note that 
this evidence is confined to the earlier stages of stellar 
evolution. 

Finally, to check the sensitivity of our results to re- 
maining uncertainties in the exponents of the scaling of the 
WR mass-loss rates, we study two extra scaling schemes 
listed in Tabled By comparing schemes B, C and D, we 
find that the rate at which the WC/WN ratio decreases is 
primarily determined by the exponent of the M — Z de- 
pendence during the WN phase. Although the error bars 
of the observed data-points are relatively large, agreement 
with the observations favours the Z-dependent schemes B 
and C. We note that it is only the quantity of the final WR 
mass that depends significantly on the ovlM — Z exponent 
during the WC phase. 



7. Summary, discussion & conclusions 

We have investigated how the predicted WR mass loss 
scaling with initial metallicity affects the output of evo- 
lutionary models. The inclusion of this scaling introduces 
important differences in WR masses and lifetimes, and 
produces a better agreement with observations, in partic- 
ular the WC/WN ratio. Of course the physical effects of 
rotation and binarity will be of additional importance. 

The primary effects of rotation will be to decrease the 
minimum initial mass for WR stars and to boost the WR 
lifetimes, thereby increasing the total number of WR stars 
and the number of typ e Ibc SNe progenitors. R otation 
has been investigated bv lMevnet fc MaedeJ l)2005l) . which 
indicated that initially rapidly rotating stars have a lower 
WC/WN ratio than non-rotating stars. In other words, if 
rotation were included in our models the WC/WN ratios 
would decrease relative to the values presented here. 

Binary stars will have similar effects to rotation in that 
they will lower the minimum initial mass for WR stars 
by providing extra opportunities for mass-loss via interac- 
tions. However, during the WR phase there will be little 
effect on the evolution, as WR stars have a smaller ra- 
dius than during previous evolutionary phases, and the 
radius will commonly be smaller than the binary orbit. 
Therefore, we expect that more WR stars will exist in re- 
lation to other stellar types (such as red supergiants) , but 
the WC/WN ratio will only be slightly affected. 
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So far w e have not men t ioned the importance of mag- 
netic fields. IPetrovic et al.l l)2005(l model rotation and bi- 
nary stars both with and without magnetic fields, and find 
significant differences in the behaviour of the stars. The 
most important effect of magnetic fields is to force a star 
to rotate as a solid body. This reduces the effect of rotation 
on the later stages of evolution as the star loses more angu- 
lar momentum than when magnetic fields are not included 
and angular mo mentum is reta ined in the star's core. The 
study by Petro vic et alJ l)2005j) comprises only a few stel- 
lar models concentrating on the question whether the WR 
stars can be the progenitors of long-duration GRBs. To 
obtain a more complete picture of the WC/WN ratio, the 
additional physical processes must be incorporated, how- 
ever this would require a large number of stellar models, 
and is beyond the scope of this article. 

In closing, although both rotation and binarity have 
important implications for WR evolution, these are an- 
ticipated of only secondary relevance to ensuring the 
mass-loss metallicity dependence, as far as WR lifetimes, 
WC/WN ratios, and final masses of compact remnants are 
concerned. 
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